Rheological Chaos in a Scalar Shear-Thickening Model 
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We study a simple scalar constitutive equation for a shear-thickening material at zero Reynolds 
number, in which the shear stress a is driven at a constant shear rate 7 and relaxes by two parallel 
decay processes: a nonlinear decay at a nonmonotonic rate R{p\) and a linear decay at rate Ao"2- 
Here cr\,i{$) = i~i~2 / &(t') exp[— (t — t')/n,2] At' are two retarded stresses. For suitable parameters, 
the steady state flow curve is monotonic but unstable; this arises when T2 > n and > R'(<j) > — A 
so that monotonicity is restored only through the strongly retarded term (which might model a slow 
evolution of material structure under stress). Within the unstable region we find a period-doubling 
sequence leading to chaos. Instability, but not chaos, persists even for the case ti — > 0. A similar 
generic mechanism might also arise in shear thinning systems and in some banded flows. 
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Rheochaos can be defined as the occurence of macro- 
scopic chaos [Q in a viscoelastic material at negligible 
Reynolds number. With the neglect of inertia that this 
implies, the nonlinearity must come not from the advec- 
tion of momentum (as in Navier Stokes turbulence) but 
from the constitutive behaviour of the material, which 
may include strong memory effects. Likewise for the 
chaos to be macroscopically observable (for example in 
time series data on the stress measured at fixed strain 
rate, or vice versa, in a bulk sample) a mechanism must 
be present that goes beyond the microscale chaos known 
to be present in, e.g., colloidal Stokes flow j|. 

Strong candidates for rheochaos include micellar mate- 
rials H , dense lamellar phases j| , and also dense suspen- 
sions where erratic stress response at fixed strain rate (or 
vice versa) is widespread but poorly documented (see e.g. 

It is not yet clear whether spatial as well as temporal 
inhomogeneity is present for all instances of rheochaos, 
and if so to what extent. This could range from a shear- 
banded flow in which the interface between bands of fast 
and slow flowing material is unsteady in time (as sus- 
pected in micelles [||6| ) through to fully developed 'elas- 
tic turbulence' as recently reported in polymer solutions 
near the overlap threshold 0. Spatial inhomogeneities 
are also known to occur in shear-thickening colloid so- 
lutions |]J^]. However, the closely related phenomenon 
of director chaos in sheared nematics has been studied 
theoretically and does not seem to require spatial inho- 
mogeneity H . In the present state of understanding, a 
theoretical search for temporal rheochaos in spatially ho- 
mogenous models remains justified. 

Recent work by the authors has studied the onset of 
temporal instability in spatially homogeneous mesoscopic 
models of shear-thickening type |l(J] . One interesting pre- 
diction was that such instability could arise in a system 
where the steady state flow curve, a(j) is monotonic p0[ . 
This contrasts with the conventional instability to spa- 
tial inhomogeneity in the form of shear bands: this is 



always associated with regions of negative slope on the 
flow curve [ [TT| ~ 13| . The mesoscopic models of 1 10 are not 
fully tensorial but work with a single (spatially uniform) 
component of each of the stress and strain rate tensors 
(cr and 7); nonetheless they contain an infinite number of 
degrees of freedom, corresponding to the distribution of 
local strain variables for different mesoscopic elements. 
This makes them complex to analyse. 

In this paper we propose closely related but much sim- 
pler models in which there is only one degree of free- 
dom (the shear stress a) whose time evolution at con- 
stant strain rate 7 is governed by a simple constitutive 
equation with retarded and nonlinear features. The sim- 
plest such model combines a nonlinear instantaneous re- 
laxation rate for stress (chosen nonmonotonic) with a 
linear but retarded relaxation. For a single exponential 
retardation kernel, its dynamics can be completely un- 
derstood: it shows spontaneous oscillation in a region of 
the flow curve with positive slope, but no chaos. This is 
qualitatively very like the mesoscopic model of Ref. |l(J 
(although that model exhibits oscillations at a constant 
imposed stress rather than strain rate). In particular, 
the instability is associated with a negative slope on the 
'bare' flow curve (before the retarded term is added). A 
second, similar model, in which the nonlinear relaxation 
is itself delayed, shows chaos. 

We first examine the simplest model alluded to above. 
This is defined by the equation 



&{t) = 7 - R(a) - \<t 2 



(1) 



where 02 (i) = /__, Ma(£ — t')a(t') dt' is a retarded stress 
and M2(t) is a memory kernel whose integral is unity. 
The first term on the right side of this equation means 
that, in the absence of relaxation, stress increases linearly 
with straining (the elastic constant is set to unity) - a 
Hookean solid. The second term describes instantaneous 
decay of stress at rate R(a), for example through 'hops' 
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or plastic rearrangement of mesoscopic elements (return- 
ing these to an unstrained state) with jump rate Rj a. 
Unlike in the mesoscopic models of |I3], no attempt is 
made to track the dynamics of individual elements. The 
third term is also a decay term, but describes retarded 
relaxation. This could represent 'delayed jumps' which, 
perhaps because they involve a cooperative motion of 
many elements, take a distribution of finite times to ac- 
complish (governed by the kernel M 2 ). More generally, a 
retarded term could represent some other slow structural 
reorganization of the material in response to stress. 

For example, one could have a model of instantaneous 
jumps but with a 'fluidity' or jump rate that itself adapts 
slowly to stress [Q . In this context it might be more nat- 
ural to have a nonlinear retarded term such as 

&=j- R(a) - \<7 2 <t (2) 

However, this gives qualitatively the same instability as 
described below for Eq.(Q) |l5[ ]; we retain the linear ver- 
sion, for simplicity, below. 

Solving Eq.(|l]) in steady state gives immediately the 
flow curve, or rather its inverse: 

7 = R(a) + Act (3) 

The interesting case is when R(a) is nonmonotonic but 
R(u) + Act is monotonic. Then the flow curve is mono- 
tonic, but only because of the retarded contribution to 
the jump rate. One might suspect that a sufficiently slug- 
gish retarded contribution might fail to correct the un- 
derlying instability in the region where R'(cr) is negative: 
over short timescales the system appears to unstable with 
respect to shear banding but at long timescales it is not. 
Here, the timescales are measured relative to the strain 
rate at which R(a) in (^|) first becomes non-monotonic; 
we choose units so that this is O(l). 

We analyse the case of a single exponential kernel, 
M2 = t^ 1 cxp[— (t — t')/T2\- As is easily checked, for this 
kernel Eq.(Q) can be replaced by a differential equation 
of second order. Differentiating Eq.(|l|) with respect to t, 
and noting that a 2 = (ct — (J2)/t2, we obtain immediately 

dV 

a = -^-i{a)a (4) 

which effectively describes a particle of unit mass in a 
1-D potential V with damping constant £. Here 

T 2 V(cr)= [ R(crf)daf + Xcr 2 /2-ia (5) 
Jo 

£(ct) = R'(a) + l/r 2 (6) 

As 7 is varied, the steady state flow curve cr( A f), as given 
by Eq.(||), is recovered as the solution of V'(a) — 0. Sta- 
bility of the steady state solution requires that two fur- 
ther conditions are satisfied. The first is V"(a) > (so 
that the effective potential has a minimum not a max- 
imum). This is equivalent to da/dj > which is the 



usual criterion to avoid shear banding. However, stabil- 
ity also requires that £(ct) is positive at the minimum of 
V. When R'(a) in Eq.(||) is negative, this is only satis- 
fied if the retardation time r 2 is sufficiently short. When 
not satisfied, one has antidamping at the minimum of 
V so that small velocity fluctuations are amplified; this 
is reminiscent of a van der Pol oscillator pq] . Velocity 
fluctuations will grow until a limit cycle is reached in 
which the postive damping at large amplitudes balances 
the antidamping near the minimum. 

An example of the 'bare' flow curve, the final flow 
curve, and the region of instability is shown in Fig. 1(a). 
Fig. 1 (b) shows typical time series of the stress just in- 
side, and well within, the unstable region. The limits of 
this region, ct^, arc Hopf bifurcation points where there 
is onset of finite frequency sinusoidal oscillations with an 
amplitude varying as I7 — A f c \ 1 ^ 2 - 

Our choice of an exponential kernel is nongeneric: 
most integral kernels are not equivalent to any finite or- 
der differential equation |l6| ]. However, the above ar- 
gument gives a generic mechanism of instability. If the 
flow curve is monotonic only because of a retarded term 
(—A < -R'(ct) < 0), then temporal instability survives 
if the retardation time is too long. Its presence does 
not depend on details of the kernel, but what it leads 
to might do so: in particular, chaos is impossible in a 
second-order system [|l6| such as Eq.(|j). However, our 
finding of spontaneous oscillation but not chaos appears 
to be structurally stable: we were unable to find chaos 
with M 2 taken as the sum of two exponentials (which 
gives a third order dynamical system for which chaos is 
allowed). 

In that case, what does need to be added to the model 
of Eq.(|]) to give temporal chaos rather than just spon- 
taneous oscillation? So far, the simplest variant we have 
found that definitely shows chaos is the following: 

&(t)=<y-R(<T 1 )-\a 2 (7) 

where the stress in the nonlinear term, cti, is now also 
retarded. The steady state flow curve is the same as that 
for (jl]). For simplicity, we choose a single exponential 
kernel here too: ct \ (£) = f Q a{t')r^ 1 exp[— (t — i')/ r i] dt' . 
To maintain continuity of interpretation with the sim- 
pler version of the model, we choose Ti ~ 1 -C T2. We 
study the situation where the monotonicity of the flow 
curve (still given by Eq.(||)) is restored only via the more 
retarded of the two relaxation terms. While there is no 
longer a simple interpretation in terms of an effective po- 
tential or a damping function, the generic instability of 
the previous model remains. But now, within the un- 
stable region, we find a period doubling cascade leading 
to chaos. Figure 2(a) shows, for a specified set of model 
parameters, the period and Lyapunov exponents Ai > 
A2 > A3 as a function of 7 (Ai > means that nearby 
trajectories exponentially separate (l^]); Fig. 2(b) shows 
a series of period-doubling orbits in the (ct 2 , ct) plane and 
Fig. 2(c) shows the strange attractor in (cti, ct 2 , cr)-space. 
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Its Lyapunov dimension -Di yap = 2 + Ai/ 1 A3 varies with 
the parameters but is slightly greater than 2, typically 

2.0 < Ayap < 2.1. 

Physically it is not clear to us yet why retardation of 
the nonlinear term (as well as the linear one) seems nec- 
essary to get chaos out of Eq. ((?]); presumably however, 
this adds something which is missing even from the meso- 
scopic model of Ref. (where chaos remained absent 
despite the infinite order of the system). Attempts to 
associate the retarded stresses in this model with, say, 
higher moments of the distribution of local strains in the 
model of Ref. (where the first moment is the instan- 
taneous stress) have so far proved unconvincing. A more 
detailed study is left for future work. 

We conclude with a broader discussion. The key idea 
is that of a flow curve (for spatially homogeneous states) 
whose monotonicity is rescued only by a retarded con- 
tribution; if too retarded, this does not restore temporal 
stability because the system continues to amplify pertur- 
bations over short timescales. Although the equations 
involved will look rather different, very similar physics 
could arise in materials of shear-thinning type where 
shear banding is present |6j|l^,[l8| or narrowly avoided 
p9| . It might be very interesting to look more closely in 
shear-thinning micellar systems where, by varying den- 
sity and temperature, one can arrange a material whose 
flow curve is only just monotonic p8| . Similar studies 
in colloids close to the transition from continuous to dis- 
continuous shear thickening || would also be valuable 
although this field is a lot less developed experimentally. 

Quite similar equations, but with different variables 
and interpretation, might describe a pre-existing shear 
banded flow, whose stability remains unclear in many 
cases J2(J. The simplest scenario would ascribe a single 
coordinate to describe the bands (e.g. the position of the 
interface between them, assumed flat) and seek to de- 
velop equations for its time evolution. Chaotic behavior 
of such an interface, rather than of a spatially homoge- 
neous stress, might be the explanation of rheochaos seen 
in various micellar systems I3J. In the case where one of 
the bands is a static gel, empirical models such as those 
proposed in Rcfs. pll] have met with some success at ex- 
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plaining the observed (though not entirely steady |22|) 
dependence of stress on strain rate when averaged across 
such a banded flow. Such models involve equations such 
as h = f(h) — 1/a where h is the width of a shear band, / 
is a nonlinear term arising from the difference in concen- 
tration in the two bands and a is the stress [^T|. Under 
controlled strain rate conditions (say) 1/a is linear in h 
and the equation is not dissimilar to Eq.(|l]) without re- 
tardation. If a slow process can be identified (possibly 
concentration equilibration), then a retarded version of 
this type of equation could share the generic instability 
of the models discussed above. 
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FIG. 1. CAPTION: (a) The bare flow curve 7 = R(a) 
(light line) and the final flow curve (heavy line) in the 
model of Eq.(^j). Parameters are A = 20, T2 = 10 and 
R(a) = 0.6(J J - 3.3ct 3 + 5(J. The region of instability 
(Tc < <y < at is shown, where o^T 0.799 and ~ 1.631, 
corresponding to ~ 18.487 and jj ~ 33.382. Note that for 
our choice of parameters, cr* almost coincide with the turning 
points of R. (b) Stress time series (same parameter values) 
at (from bottom to top) 7 = 18.49, 7 = 30 and 7 = 33.38. 
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FIG. 2. (a) Upper plot: the period of stable orbits as a 
function of strain rate 7 around the unstable region of the 
flow curve of Fig.l, for the model of Eq.(0) with ti = 0.5 and 
the other parameters as in Fig. 1. Lower plot: Lyapounov ex- 
ponents for trajectories, showing Ai > = A2 in the chaotic 
regions, (b) Orbits projected onto the (02, <r) plane for vari- 
ous 7 showing the period-doubling cascade with periods 1, 2, 
4, 8, 16 and 32. (c) The strange attractor in (<ti,(J2,(t) space 
for 7 = 20 over a time period 5 x 10 2 < t < 10 3 (arb. units). 



